Tetramerisation of a frustrated spin-1/2 chain 
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We investigate a model of a frustrated spin-1/2 Heisenberg chain coupled to adiabatic phonons 
with a general form of magnetoelastic coupling. For large enough frustration and lattice coupling 
a new tetramerised phase with three different bond lengths is found. We argue that the zig-zag 
spin-1/2 chain LiV205 might be a good candidate to observe such a phase. 
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Quasi one-dimensional (ID) quantum antiferromag- 
nets exhibit fascinating magnetic properties at low tem- 
peratures. Inorganic compounds such as CuGeOa (Rcf.|l|) 
or LiV205 (Ref.§) are almost ideal prototypes of the 
spin-1/2 frustrated chain, the so-called antiferromagnetic 
(AF) Heisenberg J1 — J2 chain (see Fig. |l|). The chem- 
istry of these compounds enables the synthesis of single 
crystals much larger than their organic analogs and con- 
sequently the achievement of new experimental studies. 
Recently, the discovery of a spin-Peierls (SP) transition 
in CuGeOa (Ref J|) has drawn both experimental and the- 
oretical interest. 

At temperatures larger than the interchain couplings 
the quasi-lD compounds CuGeOa or LiV205 are well de- 
scribed as independent AF Heisenberg chains including 
next-nearest neighbor (NNN) interactions responsible for 
frustration. The nearest neighbor (NN) Ji and NNN J2 
exchange integrals can be estimated by a fit of the mag- 
netic susceptibility, the high temperature behavior being 
governed by Ji and the position of the maximum by the 
frustration ratio J2/J1. Values such as Ji « 160 K and 
J2/J1 ~ 0.36 have been proposed for CuGeOa (RefJJ). 
On the other hand, in LiV205 the spin-1/2 V^"*" ions 
form double-chains similar to Fig. |l|(b) well separated by 
inert double-chains of V^"*" ions. Quantum chemistry cal- 
culations suggest that J2 could even be larger than Ji in 
that case |5|. 

The SP transition is an instability due to magnetoe- 
lastic effects which is characterized (below a critical tem- 
perature Tsp) by the opening of a spin gap and the ap- 
pearance of a lattice dimerisation. It was first predicted 
to occur in the non-frustrated S=l/2 chain but the 
properties of the J1 — J2 chain suggest that it is also a nat- 
ural instability in that case since the ground-state of that 
model is spontaneously dimerized for J2/J1 ^ 0.24. This 
is particularly clear at the so-called Majumdar-Ghosh 
point (MG) J2/J1 — 0.5, where the ground state 
(GS) is two-fold degenerate, corresponding to two pos- 
sible dimerisation patterns formed by a succession of dis- 
connected singlet dimers. However, when J2/ Ji becomes 
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FIG. 1: (a) Linear Ji — J2 chain; (b) Zig-zag Ji — J2 chain. 



very large, another instability could occur: The J2 chains 
are only weakly coupled, and they could undergo a SP 
transtion of their own. The interplay between both in- 
stabilities has not been considered so far. 

In this Letter, we investigate on equal footings the role 
of the frustration and of the lattice coupling. Special em- 
phasis is put on the search for new phases which would 
result from the combination of both effects. The compe- 
tition between various orderings which could eventually 
appear simultaneously can only be addressed by going 
beyond pertubative approaches. Using Exact Diagonal- 
isation techniques we report evidences for a new mixed 
phase with both dimerisation and tetramerisation ampli- 
tudes. Lastly, we discuss our results in the context of the 
quasi-lD antiferromagnets CuGeOs and LiV205. 

The Hamiltonian of a frustated spin chain on L sites 
coupled to (adiabatic) lattice displacements is written as, 

i i 

+ J2Y,i^-A2{S,+6,+,)]S,-S,+2, (1) 

i 

where Si is the distortion of the bond between site i and 
i + 1, K the spring constant and the first term corre- 
sponds to the elastic energy loss. In general, this term 
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FIG. 2: Generalised dimer susceptibilities as a function of 
chain momentum q calculated by ED on a 32 site ring. Some 
peak amplitudes for L — 20 are also shown as indicated on 
the plot, (a) NN dimer susceptibility (a = 1); (b) NNN dimer 
susceptibility {a — 2). 



might also contain cross-terms such as SiSi+i (depend- 
ing on the underlying geometry of the structure). We 
have checked that they do not affect the basis physics of 
this model so that we omit them for simplicity. Unless 
specified otherwise, Ji sets the energy scale. The spin- 
lattice couplings Aa are assumed to be dimensionless so 
that the distortions Si are given in units of the lattice 
spacing. Note also that, as can be seen from a trivial 
re-definition of the Si, the coupling strengths can be re- 
defined by the reduced variables Ai = Ai{Ji/ Ky/^ and 
A2 = A-i^JxIKY^"^ and used to investigate the phase 
diagram. However, the "physical" values for the mod- 
ulations Si depends on the A^'s and K separately. Es- 
timations of these parameters can be given on physical 
grounds as will be discussed later on. Values such that 
A2 = A\ and A2 = 2Ai are relevant for the physical 
systems we are interested in. 

The phase diagram of the frustrated chain in the ab- 
sence of lattice couplings {Ai = A2 — 0) is well known. 
The GS is uniform for small frustation (with power-law 
decay of the spin correlations) and becomes dimerized 
(with a finite spin gap) for J2/J1 larger than a critical 
value jc which has been determined with great ac- 
curacy by numerical methods, jc — 0.241167 (Ref. H). 
Interestingly enough, incommensurate spin correlations 
(away from the AF momentum q = tt) appear for 
J2/J1 > 0.5 (Ref. |l^). 

Before investigating the full Hamiltonian (|l|), it is in- 
structive to first consider "generalised" dimer suscepti- 
bilities of the form, 



Sj+aexp{iqrj))^)^, (2) 



where the expectation value (...)g is taken in the GS 



of the J1—J2 chain in the absence of lattice coupling. 
Physically, any instability towards a modulated dimer 
phase involving dimers at distance a would be signaled 
by a sharp peak of S {q, a) at a given q associated to 
the wavevector of the modulation. As seen in Fig. ^ 
(Ref. sharp peaks are indeed seen in S{q,l) and 
S{q,2) at momentum q — t: and q = 7r/2 respectively 
signaling proximity of instabilities toward the formation 
of dimerised {q = vr) and tetramerized {q = tt/2) phases 
involving NN and NNN dimers respectively. Note that 
NNN g = 7r/2 dimer correlations increase with increasing 
frustration J2 / Ji while the maximum NN dimer suscep- 
tibility occurs around J2/J1 — 0.5. Since the order pa- 
rameter {^j- Sj ■ Sj^2 exp («f J'j)) of the tetramerised 
phase is directly coupled to A2 in Hamiltonian the 
finite magnetoelastic A2 coupling is then the key feature 
of the model. 

In order to solve Hamiltonian (|^) including the mag- 
netoelastic coupling we use Lanczos diagonalisations of 
small finite rings of size L with periodic boundary con- 
ditions. An iterative procedure is used to determine the 
displacements Si by solving a set of coupled non-linear 
equations [T2| , 

KS, - JiAi{S^ ■ S^+l) (3) 
- J2A2{{S, ■ S,+2) + (5^.-1 • 5,+i)) = 0, 

where (...) is the expectation value in the GS of Hamil- 
tonian (|l|). Note that no translation symmetry is a pri- 
ori assumed in order to search for lattice modulations 
of arbitrary periodicity (compatible with system size). 
We found that, generically, the lowest energy lattice con- 
figuration is obtained for a single or a superposition of 
the following distortions; (i) a uniform (negative) com- 
ponent Si = So (which is due to the finite compressibility 
of the system |l^), (ii) a dimerisation Si = (5d(— 1)' and 
(iii) a tetramerisation Si = St cos {^i + (I)t)- Although, 
the tetramerisation could be either site-centered (with 
(j)rp = 7r/4) or bond-centered (with (j)^ — 0), only the sec- 
ond bond-centered type (i.e. a modulation of the bonds 
like St-0-{—St)~0) was found. This particular pattern 
can easily be understood in the large- J2 limit which con- 
sists of two weakly coupled J2 Heisenberg chains (see 
Fig. |l|(b)). In that limit, the magnetoelastic coupling A2 
tends to produce a dimerisation of each chain so that 
{S2p + S2P+1) tx (—1)'' which can indeed be realized by a 
tetramerisation of the chain with 0t = 0. 

The domains of stability of the various phases are 
shown in Fig. |^ for two values of the parameter ^2/^1- 
Data are shown for cyclic rings of size L = 12, 16 and 
20 so that tentative phase diagrams can be obtained (see 
Fig. ^ from a finite size scaling analysis. Various regimes 
have to be distinguished for these extrapolations. First, 
when J2/J1 < jc, the Heisenberg chain is critical and 
one expects that a dimerised GS would be stabilized 
for arbitrary magnetoelastic coupling |q|. Indeed, the 
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FIG. 3: Phase diagram in the A\ - frustration plane for dif- 
ferent system sizes L. The data points (and the dotted lines) 
correspond to the critical values of the (reduced) magnetoelas- 
tic coupling above which dimerisation and/or tetramerisation 
(as shown on the plot) appear; (a) A2I A\ = 1; (b) A2lA\ = 2. 
In the latter case, a very small region of uniform phase inside 
the dimerised phase area (omitted here for clarity) was found 
but shown to be a spurious finite size efi'ect. 



finite critical value of the coupling exhibits a clear 
Xj \Il power-law behavior with system size. In the range 
jc < J2I J\ < 0.5 the system size dependence of the crit- 
ical value of Ai becomes exponential signaling the fact 
that the infinite Heisenberg chain forms singlet dimers, 
even in the absence of the lattice. Note that at the spe- 
cial MG point J2/J1 = 0.5 the critical coupling for Ai 
vanishes for all sizes. Special care is needed to analyse 
the data for J2/J1 > 0.5: In a narrow range of J2/J1, 
0.5 < J2/J1 < jcrit,n which depends on the A2/A1 ratio 
the system only dimerises above a small critical value of 
the coupling Ai. In a range jcrit,i < J2/J1 < icntai^) 
(which extends with increasing system size), when Ai 
exceeds a higher critical value, a tetramerisation super- 
poses to the existing dimerisation. For J2/J1 > Jcrit,2(i) 
dimerisation and tetramerisation occur both for the same 
critical value of the coupling constant. However, our 
data are consistent with the fact that Jcrit,2(i) ^ 00 
when i — > cx) so that this last regime seems irrelevant. 
For 0.5 < J2/J1 < jcrit,2(-^), the critical coupling for 
dimerisation rapidly vanishes with increasingly large sys- 
tem sizes (although some increase has been observed for 
small sizes). For J2/J1 > jcrit,i (e.g., icrit.i ~ 0.85 for 
A2 — Ai), the finite size dependence of the critical cou- 
pling for tetramerisation is consistent with a rapid ex- 
ponential behavior converging to a finite value. An es- 
timation of the infinite size phase diagram is then pos- 
sible with reasonable accuracy as shown in Fig. ^. In 
summary, our calculations predict that the J1 — J2 chain 
is always dimerised once it is coupled to the adiabatic 
lattice. Tetramerisation, however, occurs only for large 
enough frustration and when the lattice couplings Ai and 
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FIG. 4: Tentative phase diagrams in the Ai - frustration 
plane obtained from a finite size scaling analysis of the data of 
Fig. ^ The dimerized and tetramerized phases are labelled by 
D and T respectively. Thick dashed (full) lines correspond to 
first order-like (second order) transition lines, (a) A2/A1 — 1; 
(b) A2/A1 = 2. 



A2 exceeds some critical values which vanish when J2/ Ji 
increases to infinity. Interestingly enough, our data also 
suggest that the nature of the D ^ T transition might 
change when frustration increases, from first order (with 
discontinuities in the dimerisation and tetramerisation 
amplitudes) to a continuous second order-like line at large 
J2. Note that, even for large J2/J1 (limit of weakly cou- 
pled chains), one still expects finite critical lattice cou- 
plings since, in the absence of the lattice, an arbitray 
small perturbation Ji introduces an exponentially small 
gap and a finite spin correlation length. 

We finish by discussing some applications of the 
present study to real materials. The magnetoelastic cou- 
plings are generically due to strong dependence of the 
exchange integrals with respect to distance, typically 
Ja(r) cx r~"", with an exponent aa in the range 7-15. 
A small change of length Sr of the bond connecting two 
sites at distance along some direction Ua (r = raUa) 
leads to a linear change of the AF coupling, 



Ja{5r) = Ja(l - —Ua 



Sr) 



(4) 



In the case of the linear chain of Fig. |i|(a) where the dis- 
placements occurs along the chain direction, Eq. (^ pre- 
dicts A2 — Ai/2 assuming the same values of aa for the 
two chemical bonds. In the case of CuGeOa, the superex- 
change path giving rise to J2 involves more intermediate 
states (in particular Germanium orbitals) so that one ex- 
pects a2 > oil and A2 ~ Ai seems more physical in that 
case. According to the phase diagram of Fig. |^(a), condi- 
tions for a small tetramerisation seem clearly not realised 
in CuGeOa. Indeed, for a frustration J2/J1 ~ 0.4 and 
a small physical value for the dimensionless coupling Ai , 
the A2 coupling becomes irrelevant (apart from produc- 
ing a tiny overall contraction of the lattice) and we expect 
a simple dimerisation 60 oc AiJi/K . Assuming Ai ^ 10, 
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FIG. 5: Amplitudes of the relative bond length change \5o\ 
(open circles), dimerisation Sd (diamonds) and tetramerisa- 
tion 5t (stars) (in units of the lattice spacing) versus frus- 
tration calculated on L = 20 site rings for Ai = 5 and 
A2 — 10. Only non-zero amplitudes are shown; (a) Ai = 1; 
(b) Ai = 1/3: data for L = 16 (thick dot-dashed lines) are 
also shown in that case to indicate finite size effects. 



To conclude, from numerical calculations we have ob- 
tained the generic properties of the frustrated spin-1/2 
chain coupled to adiabatic phonons. Our results are 
confronted to experimental systems such as the (quasi) 
linear CuGeOs chain and the zig-zag LiV205 chain. 
While the observed small dimerisation of the SP phase of 
CuGeOa can easily be accounted for within such a simple 
model, we argue that LiV205 is a good candidate to ob- 
serve a novel tetramerised phase with a doubling of the 
zig-zag chain periodicity (leading to new superstructure 
peaks in diffraction experiments) and two (slightly) non- 
equivalent V''+ sites, a feature which could be observed 
in NMR-experiments. 
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Ji ~ 100 Kelvins and K ^ 10 eV a dimerisation around 
0.1-0.3% of the lattice spacing is expected in agreement 
with X-ray diffraction experiments p^ . 

We now turn to the case of the LiV2 05 compound. If 
one assumes that atomic displacements in Fig. |l|(b) would 
occur along the zig-zag chain direction and that the expo- 
nents aa are identical for the two bonds, Eq. (jj) implies 
that A2 = 2Ai. As shown in Figs. ||(a-b), the stability 
of the tetramerised phase increases with increasing ratio 
A2/A1. Hence, due to larger A2/A1 and J2/J1 ratios, 
LiV205 seems, contrary to CuGeOa, an interesting can- 
didate for the new tetramerised phase. Typical exponents 
like ai = a2 = 10 gives Ai = 5 and A2 = 10. Assum- 
ing a physical value Ji = 400 Kelvins for the exchange 
constant and choosing K ~ IcV {— 10 000 Kelvins) 
and K ^ 9eV 90 000 Kelvins), we get ii = 1 and 
Ai = 1/3 respectively. The magnitudes of the dimerisa- 
tion and tetramerisation for these parameters are shown 
in Fig. H as a function of the frustration J2 / Ji . For rather 
large magnetoelastic couplings such as the one used in 
Fig. ||(a) where finite size effects are negligible we ob- 
serve, for increasing magnetic frustration, a transition 
from a purely dimerized phase to a new phase with a 
dominant tetramerisation and a small dimerisation com- 
ponent. For realistic couplings, let's say K/Ji > 200, as 
seen m Fig. |(b), finite size effects become large. Never- 
theless we expect a behavior similar to that of Fig. ||(a) 
although with much smaller lattice displacements. Typ- 
ically, while 6t might be of the order of a percent of the 
lattice spacing, Sd is expected to remain much smaller. 
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